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Simulation and Analysis Chain for Acoustic Ultra-high 
Energy Neutrino Detectors in Water 
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Abstract. Acoustic neutrino detection is a promising approach! for large-scale ultra-high energy neutrino detectors in water. In 
this article, a Monte Carlo simulation chain for acoustic neutrino detection devices in water will be presented. The simulation 
chain covers the generation of the acoustic pulse produced by a neutrino interaction and its propagation to the sensors within 
the detector. Currently, ambient and transient noise models for the Mediterranean Sea and simulations of the data acquisition 
hardware, equivalent to the one used in ANTARES/AMADEUS, are implemented. A pre-selection scheme for neutrino-like 
signals based on matched filtering is employed, as it is used for on-line filtering. To simulate the whole processing chain for 
experimental data, signal classification and acoustic source reconstruction algorithms are integrated in an analysis chain. An 
overview of design and capabilities of the simulation and analysis chain will be presented and preliminary studies will be 
discussed. 
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PACS: 95.55.Vj, 95.85.Ry, 43.30.Zk 



1. INTRODUCTION 

Acoustic neutrino detection uses the effect that ultra- 
high energy (UHE, E^ > lO^GeV) neutrinos can pro- 
duce a detectable acoustic pulse according to the thermo- 
acoustic model [1]. This model describes the genera- 
tion of an acoustic pulse due to the local heating of the 
medium by a particle shower, which is caused by a UHE 
neutrino interaction. The fast energy deposition by the 
cascade in the medium leads to a bipolar pulse. Due to 
coherent emission over the cylindrical geometry of the 
cascade, the wave propagates through the medium in a 
disk-like shape perpendicular to the main axis of the cas- 
cade. Given the expected low flux of neutrinos with en- 
ergies in excess of lO^GeV, a potential acoustic neu- 
trino telescope must have large dimensions of presum- 
ably > lOOkm^. The attenuation length of sound in wa- 
ter is of the order of 1 km for the peak spectral density of 
a UHE neutrino induced signal of around 10 kHz, allow- 
ing for a sparse instrumentation. In this article, a simu- 
lation and analysis chain for an acoustic neutrino detec- 
tor in water is presented. This includes the reproduction 
of the acoustic pulse generation, the detector properties, 
and the deep-sea acoustic environment. To simulate and 
evaluate the complete processing chain for experimental 
data, strategies for the position reconstruction and classi- 
fication of acoustic signals were integrated. Preliminary 
results for an effective volume and a transient-free limit 
setting potential of the AMADEUS [2] acoustic detection 
system will be discussed. 



2. SIMULATION CHAIN 

The simulation chain [3] consists of the following mod- 
ules, which build upon each other to create a simulated 
event: 

• An interaction vertex is located at a random posi- 
tion in a given volume around the detector and the 
energy and direction of the shower are set randomly 
within definable ranges. 

• The formation of the shower and the resulting 
acoustic signal, as generated by a UHE neutrino in- 
teraction, are simulated. 

• The acoustic environment of the deep sea is re- 
produced including both the ambient and transient 
noise conditions. 

• The data acquisition (DAQ) hardware is simulated 
including the system response and inherent noise of 
the sensors and read-out electronics. 

• The pre-selection scheme used on-line in the detec- 
tor for data reduction is applied. 

The Monte Carlo (MC) shower is produced from a 
parametrisation, which is based on work by the ACoRNE 
collaboration [4, 5]. This parametrisation, which is valid 
up to a total shower energy of lO^^GeV, describes the 
distribution of the deposited energy in the surrounding 
water After the cascade has been simulated, the acoustic 
pulse and its propagation to the sensors within the detec- 
tor are calculated. The deposited energy of the shower 
produces a local heating of the medium. With respect 



to hydrodynamical time scales, the energy deposition at 
time to is instantaneous. The energy deposition e(r,f) 
can be factorized into a spatial and temporal part, where 
the latter is approximated by using a Heaviside function. 
Assuming a total energy deposition £ in a cylindrical 
volume around the shower, the spatial part e(r) can be 
expressed for the longitudinal and radial positions z and 
r in the shower as: 



e(r) 



1 1 d^E 
E 2nr drdz 



(1) 



An expression for the pressure p can be derived, follow- 
ing [6]: 



P{r,t) = 



Ea 

AnC„ 



dV ,,, d^, R. 
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(2) 



where R = |r — r'| can be approximated by the distance 
between the shower maximum and the sensor, a is the 
thermal expansion coefficient, Cp is the specific heat ca- 
pacity at constant pressure and Cj the speed of sound 
in water From the parametrisation of the spatial energy 
density distribution, which is deposited by the cascade, 
the pressure at a sensor at a distance R from the shower 
can be calculated numerically [3, 5]. The pressure signal 
is then convoluted with the frequency dependent sound 
attenuation in sea water. This attenuation is based on a 
model by Ainslie and McColm [7], extended to a com- 
plex representation. This complete procedure calculates 
a UHE neutrino induced bipolar pulse at the sensor in- 
cluding the disk-like propagation pattern. In Fig. 1, an 
example of a simulated pressure pulse is shown. 
The background for acoustic neutrino detection in the 
deep sea consists of two different types of noise: tran- 
sient and ambient noise. Transient noise signals have 
short duration and an amplitude that exceeds the ambi- 
ent noise level. These signals can mimic bipolar pulses 
from neutrino interactions. Sources of transient signals 
can be marine mammals and anthropogenic sources such 
as shipping traffic. The ambient noise is mainly caused 
by agitation of the sea surface [8], i.e. by wind, break- 
ing waves, spray, and cavitations. Thus it is correlated 
to the weather conditions, mainly to the wind speed [9]. 
The model used for the simulation of the ambient noise 
is based on the so-called Knudsen spectra [10], which are 
adapted to the deep sea by applying attenuation effects. 
The wind speed distribution included in the simulation 
was derived from measurements at several weather sta- 
tions near the coast of Marseilles, France. As reproduced 
by the simulation, the mean noise level (cTnoise) is about 
25 mPa for the frequency range from 1 — 100 kHz and for 
95 % of time the noise level is smaller than 2((7noise) [1 !]• 
The simulation of the DAQ hardware comprises two 
parts: the simulation of the sensors and of the read-out 
electronics. The design of the DAQ hardware is inspired 
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FIGURE 1. Simulated pressure signal for a 10 GeV 
shower as recorded at the position of the sensor. The connection 
line between the sensor and the shower maximum is perpen- 
dicular to the main axis of the shower; the distance between 
shower and sensor is 1000 m. 



by the AMADEUS [2] project. This includes acoustic 
sensors using the piezo-electric effect (hydrophones) and 
read-out electronics to amplify and digitise the signal. 
The inherent noise and the system transfer function for 
both the sensors and the electronics have been measured 
in the laboratory. Sensors normally show a directional 
dependency of their sensitivity, therefor signal and am- 
bient noise have to be treated separately. In this case, the 
incident direction of the noise is the sea surface above 
the detector Signal and ambient noise are then superim- 
posed. The inherent noise of the sensor is added. The re- 
sulting waveform is convoluted with the system transfer 
function of the read-out electronics and the correspond- 
ing inherent noise is added. 

The output is directed to a simulation of an on-line filter 
system [12], which, for real data, is used to reduce the 
amount of data to store and to pre-select waveform sam- 
ples for further off-line analysis. The filter is based on a 
matched filtering technique using a pre-defined bipolar 
pulse as reference to select signals with bipolar shape. 
In addition, a coincidence test between the sensors of a 
cluster can be performed. 



3. ANALYSIS CHAIN 

The simulation chain is used to study signal classifica- 
tion [13] and position reconstruction algorithms [14]. For 
the presented studies, a detector geometry identical to 
the configuration of the AMADEUS detector [2] is used, 
assuming six clusters of six sensors each at fixed posi- 
tions and with fixed orientation. Different sensor spac- 
ings are implemented ranging from 1 m within clusters 
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FIGURE 2. Sketch of the geometry of the AMADEUS de- 
tector. Sensor clusters are indicated by oval shapes and are 
arranged along two vertical structures. Adapted from [2]. 
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FIGURE 3. Angular resolution of the direction reconstruc- 
tion shown for the zenith and azimuth angle. The mean of the 
distribution is around zero and the sigma is about 0.7 ° for both 
zenith and azimuth angle. The acoustic sources were generated 
in a cube of 5 x 5 x 2.5 km^ around the detector centre. 



up to 350 m between clusters (cf. Fig. 2). The precise 
reconstruction of the arrival time of the signal is cru- 
cial for the direction and position reconstruction of the 
acoustic source. The arrival time is determined by cross- 
correlation with a pre-defined bipolar pulse. This proce- 
dure achieves a precision of about I jj.s. The direction re- 
construction is based on a least square fit of the measured 
arrival times at a given sensor cluster: 



min(2^(rmeasured, — f expected, ( ^ : ) ) 



(3) 



where / G 1 . .A^ is the /-th sensor of a cluster of A^ sen- 
sors, t is the arrival time, and 9 and (p are the zenith 
and azimuth angle, respectively. The acoustic sources for 
this analysis were generated in a cube of 5 x 5 x 2.5 km^ 
around the detector centre. The angular resolution ob- 
tained with the direction reconstruction algorithm is cen- 
tred around zero with a sigma of about 0.7 ° for both 
zenith and azimuth angle as shown in Fig. 3. 
The position of the acoustic source is reconstructed us- 
ing a ray tracing technique. If the directions were recon- 
structed for at least two of the sensor clusters, the in- 
tersection point (or its best approximation) of the rays 
starting from the sensor clusters and pointing into the re- 
constructed direction is searched for In Fig. 4, the distri- 
butions of the deviation of the x, y, and z-coordinate of 
the reconstructed position from the true position of the 
vertex, as used for the simulation, is shown. The half- 
width at half-maximum (HWHM) of the distributions is 
better than 15 m for each coordinate. 
The classification strategy [13] stems from machine 
learning algorithms trained and tested with data from 
the simulation. Random Forest and Boosted Trees algo- 
rithms have achieved the best results for individual sen- 
sors and clusters of sensors. For individual sensors, the 
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FIGURE 4. The distributions of the deviation of the x, y, and 
z-coordinate of the reconstructed position from the true position 
of the vertex. The acoustic sources were generated in a cube of 
5 X 5 X 2.5 km around the detector centre. 



classification error is of the order of 10% for a well 
trained model. The combined results of the individual 
sensors in a cluster are used as new input for training. 
This method obtains a classification error below 2 %. 



4. EFFECTIVE VOLUME 

The simulation and analysis chain described above was 
used to simulate the data required for this study. An 
effective volume Veff for the AMADEUS detector can be 
defined as: 



Veff(£v) 
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FIGURE 5. The effective volume of the AMADEUS detec- 
tor as a function of the logarithmic neutrino energy for the three 
different Levels, as described in the text. Also shown are the 
random coincidence rates for the Levels 1 and 2/3. 



where A^gen is the number of generated neutrino interac- 
tions in a volume Vgen and p{Ev,v,ep) is the probabil- 
ity that the neutrino reaches the interaction vertex set in 
the simulation. 5sei G {0,1} accounts for the fact that the 
probability only contributes to the effective volumeif the 
pressure pulse corresponding to the neutrino interaction 
was selected by the on-line filter within a time window of 
128 jUs around the expected arrival time. The probability 
that the neutrino reaches the vertex is given by: 



p(£'v,r,ep)=e-''wE(r,ep)A„ 



(5) 



where r is the position of the interaction vertex, Cp is 
the unit vector of the direction of the flight trajectory. 
The mean free path Awater(£v) of the neutrino in water 
is anti-proportional to the neutrino's total cross section. 
The total cross section as a function of the energy Ey was 
parameterised using values from [15]. The distance c/vve 
is the water equivalent of the distance traveled through 
matter of varying density encountered by the neutrino 
along its flight path. For the determination of the den- 
sity distribution along the flight path, the PREM' [16] 
was used to model the earth's density profile. In addi- 
tion, it is assumed that the earth is covered by water of 
2.5 km depth and the detector is placed on the sea floor 
For the calculation of the effective volume, 10^ neutrinos 
with energies uniformly distributed between lO^GeV 
and lO'^ GeV, which is the highest energy, for which the 
simulation is valid, were simulated. The uniform energy 
distribution was chosen to assure a sufficient number of 
events over the entire energy range. The interaction ver- 
texes of these neutrinos were chosen in a cylindrical vol- 
ume of 1200km^ around the AMADEUS detector The 
heading of the flight path was ranging from — 360° in 



azimuth and from — 100° in zenith^. Neutrinos enter- 
ing the generation volume from below the horizon will 
traverse an increasing amount of matter. For a zenith an- 
gle greater than 100 °, the probability of a neutrino in the 
energy range under consideration to reach the interaction 
vertex is practically zero. To determine random coinci- 
dences formed by the ambient noise, a separate set of 
simulated data was created not containing any signals. 
The model used to simulate the ambient noise was de- 
scribed in Sec. 2. The effective volume has been calcu- 
lated for three different "Levels" describing increasingly 
realistic conditions and reconstruction requirements: 

Level 1: The ambient noise is assumed to be minimal, 
matching sea state at all times, and the coinci- 
dence requirement for the filter simulation is that 
at least two sensors on one storey need to respond. 

Level 2: The complete ambient noise model and the 
standard on-line filter of AMADEUS are used, re- 
quiring at least four sensors on two storeys each. 

Level 3: In addition to Level 2, the reconstruction of 
the acoustic source position is required to have a 
deviation from the simulated position of less than 
100 m. 

The results of this study are shown in Fig. 5 for the 
three Levels. For Levels 1, 2, and 3, the effective volume 
is statistically incompatible with random coincidences 
above 1.8 • 10^" GeV, 1.8-10" GeV, and 3.2 • lO'^ GeV, 
respectively. The requirements of Level 1 are minimal, 
so this can be seen as an idealised detection threshold 
of the AMADEUS detector. The effective volume for 
Level 1 exceeds 2km^ at lO'^GeV, while for Level 2, 
it is about 0.1 km^ and, for Level 3, below 0.01 km^ 
at this energy. This result for Level 3 shows that the 
determination of the position of the interaction vertex 
has sizeable uncertainties for a small detector with an 
essentially two-dimensional configuration such as the 
AMADEUS system. 



5. TRANSIENT-FREE LIMIT SETTING 
POTENTIAL 

Following the approach derived in [17], a flux model 
independent limit estimate assuming complete transient 
noise suppression can be calculated as: 

A'90%CL 



4>, 



90%CL ■ 



(6) 



Preliminary earth reference model 



Q.TEy[V,s{Ey)/X{Ey)Y 

where N()q%cl — 2.44 for the 90% confidence level 
(CL) assuming no background and no true signal ob- 
served [18], ii is the solid angle, T is the integrated mea- 
surement time, Veff is the effective volume and A is the 



A zenith angle of ° corresponds to a neutrino coming from above. 
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FIGURE 6. The transient-free limit setting potential of the 
AMADEUS detector is shown together with the theoretical 
neutrino flux prediction [19] (pattern). For detail description, 
see the text. 



mean free path of the neutrino as described in the pre- 
vious section. For the limit estimate of the AMADEUS 
detector to UHE neutrinos, the effective volume V^s for 
Level 2 and Level 3 were used, as described in the previ- 
ous section. The assumed integrated measurement time 
used for this calculation is one year and the energy range 
and solid angle are the same as used to calculate the 
effective volume. In Fig. 6, the transient-free limit esti- 
mates of the AMADEUS detector are shown together 
with predictions for the cosmogenic neutrino flux [19]. 
The limit estimates derived show the potential of the 
acoustic UHE neutrino detection technique. Clearly, the 
limit setting potential of the detector is limited by its 
small effective volume. Complete transient noise sup- 
pression is challenging and relies on a good position re- 
construction to define a sensible fiducial volume. As the 
results for Level 3 demonstrate, the size and geometry of 
the AMADEUS detector are not favourable for this task. 



6. CONCLUSION AND OUTLOOK 

As shown, the simulation chain is capable of reproduc- 
ing the aspects necessary for acoustic neutrino detec- 
tion - from the generation of the bipolar acoustic sig- 
nal to the simulation of different detector geometries 
and components. Also the deep-sea environment with 
its variable and diverse noise conditions is well repli- 
cated. Furthermore, the analysis chain for experimental 



data, employing position reconstruction and classifica- 
tion algorithms, was verified using simulations. Based 
on this software framework, a calculation of the effec- 
tive volume and the transient-free limit setting potential 
of the AMADEUS detector was performed. For realis- 
tic conditions, the effective volume is around 0. 1 km^ at 
lO'^GeV. For a small detector such as AMADEUS, the 
calculated transient-free limit estimate is promising and 
encourages further studies for large-scale detectors. 
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